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I.  INTRODUCTION 


September  started  badly.  A  string  of  14  successful  tests  of  a  gas 
generator  at  its  cold-temperature  limit  ended  in  the  simultaneous  failure  of 
two  units.  The  meager  data  suggested  four  plausible  causes:  moisture 
intrusion,  gas  leak,  depressurization  extinguishment,  or  marginal  design.  A 
diagnostic  test  program  was  designed  to  isolate  the  most  likely  cause  among 
moisture,  leaks,  and  ignition  overdrive.  Elliott  at  the  Jet  Propulsion 
Laboratory  started  building  a  mathematical  model  of  the  generator. 

October  ended  badly.  Only  one  test  condition  (moisture)  duplicated  the 
failure.  But  all  tests  with  the  entire  generator  soaked  at  the  cold- 
temperature  limit  failed  by  dropping  below  the  minimum  acceptable  pressure. 
Although  moisture  probably  caused  the  specific  September  failures,  a  fatal 
design  flaw  had  been  exposed.  Tests  discredited  the  prime  suspect,  leakage. 

If  it  wasn't  moisture  or  leaks,  what  was  it? 

The  clearest  clue  came  from  Elliott's  model  where  the  low  pressure  could 
be  explained  by  a  steeper  than  expected  dependence  of  propellant  burning  rate 
on  pressure  at  low  pressure.  But  there  were  essentially  no  useful  data  on 
cold-temperature  burning  rate  at  that  low  pressure.  If  the  slope  were  high, 
any  drop  in  pressure  would  tend  to  continue.  Thus  when  the  cold  hardware 
cooled  the  gases,  the  pressure  drop  from  its  ignition  peak  would  not  stop  at 
the  design  limits  calculated  with  the  slope  obtained  from  higher  pressure 
burning-rate  data. 

Two  changes  were  considered:  a  new  design  or  a  higher  cold  temperature 
limit.  Schedule  demands  said  raise  the  limit,  which  had  been  -25  F,  to  where 
the  problem  disappeared.  No  one  knew  that  limit  yet.  But  when  a  conservative 
limit  of  20  F  also  failed,  design  change  awakened.  The  easiest  change  added  a 
venturi  to  isolate  the  combustion  chamber  from  pressure  loss  in  the  downstream 
tubing.  Meanwhile,  the  propellant  lot  that  had  been  used  for  two  years  had 
run  out  and  newly  made  propellant  was  being  delivered.  The  combination  of  the 
venturi  and  new  propellant  ended  the  failures. 

The  hardware  now  worked  even  at  the  lowest  temperature  but  it  was  not  yet 
clear  why.  Would  it  still  work  two  years  later?  What  were  the  variation 
limits  of  combustion?  Burning-rate  tests  validated  the  slope-break  hypothesis 
and  found  some  lot-to-lot  variation  in  the  low-pressure  slope  and  the  break 
point . 

Three  circumstances  created  the  mystery:  (1)  inadequate  knowledge  of  the 
low-pressure  burning  rate,  (2)  no  useful  simulation  of  the  generator 
operation,  and  (3)  inadequate  recognition  of  the  coupling  among  heat  transfer, 
pressure,  and  burning  rate.  What  remained  was  to  explain  the  success,  the 
failure,  and  the  probability  of  repeating  the  failure. 


II.  THE  GENERATOR  DESIGN 

The  generator  burns  an  ammonium-nitrate-oxidized  rubber  propellant  in  ;n 
end-burning  grain  to  produce  a  clean,  cool  (1300  K)  gas.  The  gas  passes 
through  a  coarse  pre-filter,  a  swirl  flow  centrifugal  filter,  a  massive  valv-; 
and  40  cm  of  metal  tubing  to  exit  a  deLaval  nozzle.  In  the  ignition  sequence, 
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a  squib-activated  igniter  spews  hot  gas  and  particles  onto  a  pellet  of  an 
ammoniura-perchlorate-based  propellant  (denoted  as  Ignition  Pellet)  which,  in 
turn,  furnishes  enough  hot  gas  to  ignite  the  main  propellant  grain  and  a 
booster  pellet  of  the  same  composition.  After  about  one  centimeter  of  burning 
of  the  main  grain,  the  propellant  geometry  becomes  a  constant-area  end-burning 
grain.  All  tubing  is  large  enough  to  keep  the  flow  velocity  below  20  m/s. 

A  venturi  of  throat  area  slightly  larger  than  the  nozzle  was  inserted 
just  upstream  of  the  valve  and  the  last  30  cm  of  tubing. 


NOZZLE 

Figure  1.  Schematic  of  Generator 


III.  MODELING  APPROACH 

Two  models  were  developed  in  this  study:  a  one-dimensional-flow  model 
and  a  lumped-parameter  model.  The  lumped-parameter  model  divides  the 
generator  into  two  chambers  separated  by  the  venturi.  It  ignores  axial 
variation  of  the  flow  except  as  the  downstream  chamber  is  separate  from  the 
upstream  chamber.  The  one-dimensional-flow  approach  addressses  the  axial 
dependence  of  the  fluid  mechanics  from  propellant  surface  to  the  exit 
nozzle.  The  one-dimensional  model  is  more  accurate,  but  expensive.  The  crude 
but  thrifty  lumped-parameter  model  allows  many  arbitrary  changes  of  input  data 
or  equations  for  sensitivity  tests.  In  retrospect  the  two  model  approach 
proved  wise. 


IV.  A  LUMPED-PARAMETER  MODEL 

The  generator  will  be  treated  as  two  chambers  separated  by  a  venturi. 
The  upstream  chamber  (denoted  1)  includes  the  propellant  combustion  chamber, 
the  filters,  and  the  tubing  up  to  the  venturi.  The  downstream  chamber 
(denoted  2)  includes  the  valve,  and  the  downstream  tubing  to  the  nozzle. 
Ordinary  differential  equations  describe  conservation  of  mass  and  energy  in 
each  chamber. 
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where  ui  ,  w  and  u  are  mass  fluxes  from  the  propellant,  through  the 

prop  ven  noz 

venturi,  and  nozzle,  respectively.  The  term  q  is  heat  loss  from  the  volume  in 
question. 


Heat  transfer  from  the  gas  to  its  bounding  surfaces  is  by  convection, 


4w  =  hc  <Tg  "  Tw>  (5) 

where  the  coefficient,  hc ,  depends  on  gas  properties  and  velocity.  Heat  loss 
to  the  surroundings  is  by  unsteady  convection  to  the  combustion-chamber  head, 
tin;  phenolic  filter,  and  the  valve.  It  is  by  free  convection  and  radiation 
from  the  tube  walls. 

The  internal  convective  coefficient  was  taken  from  a  design  analysis 
where  it  depended  only  on  mass  flow  rate,  once  all  other  conditions  were 
fixed.  Wide  variance  in  flow  velocity  from  the  reference  condition  would 
introduce  additional  error. 

h  =  110  (-— v)°‘8  [3TU/(hr-f  t2-sec)  ]  (6) 

c  .003 


where  m  is  mass  flow  rate  in  Ibm/sec.  Unsteady  conduction  may  be  approximated 
by  Go'dman's  cubic  profile  to  yield  surface  temperature. 


T  =  T  -  n  +  [ (T  -n)2  +  2nT  -  T2  ] 1/2  (7) 

wo  o  go 


where  n  =  T-  h  H/k2 
3  c 


T.R.  Goodman.  "Application  of  Integral  Methods  to  Transient  Nonlinear  Heat 
Transfer. "  Advances  in  Heat  Transfer,  Ed.  By  T.F.  Irvine  and  J.P.  Hartnett, 
Vol.  I,  Academic  Press,  1964. 
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and  a  is  thermal  diffusivity.  Radiation  is  by  the  standard, 

q  =  a  (T4  -  T4)  (9) 

r  wo 


where  the  temperature  of  the  surroundings,  T  is  the  initial  temperature. 

Burning  rate  of  the  propellant  comes  from  the  propellant-maker's  test 
data.  The  rate  is  fitted  to  the  usual  power-law  dependence  on  pressure, 


r  = 


aPn 


(10) 


This  generator  had  a  peculiarity  in  that  the  burning  rate  of  the  main-grain 
propellant  had  a  much  higher  exponent  (n)  at  low  pressure  than  at  high 
pressure  (see  Table  1).  No  problem  is  introduced  to  the  modeling  once  the 
data  are  available  and  correctly  interpreted. 

Flow  through  the  nozzle  is  assumed  choked,  quasi-steady,  and  isentropic 
and  thus  calculated  by 


0) 

noz 


1/2 


Y+l 

f  2  .  2(y-1) 

Y  +  l 


(11) 


where  Pj  and  P2  are  pressures  in  the  upstream  and  downstream  chambers, 

respectively.  is  the  cross-sectional  area  of  the  nozzle  throat.  The  gas 

n 

temperature  (T2)  used  here  was  not  the  average  chamber  temperature.  Instead, 
the  axial  temperature  drop  through  the  downstream  chamber  was  estimated,  and 
a  nozzle  inlet  temperature  calculated.  The  result  is  a  lower  temperature  and 
higher  mass  flow  rate  than  would  be  calculated  by  a  purely  lumped-parameter 
approach. 


Mass  input  to  the  upstream  chamber  from  propellant  combustion  is 


oj  =  rp  S 
prop  s 


(12) 


Flow  through  the  venturi  from  upstream  to  downstream  chamber  depends  on 
whethP’-  the  venturi  is  choked.  The  criterion  for  choking  comes  from 
isentropic  flow.  The  venturi  is  choked  (and  obeys  an  equation  similar  to  Eq . 
(11))  whenever 
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1  + 


(Ill) 
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M2] 


(-J-) 

1  — Y 
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where  M  satisfies 


Y+l 

2  (y-1 ) 


(14) 


CfM  [I  +  -  M2  ] 


is  defined  in  Eq.  (28),  and  is  the  cross-sectional  area  of  the  cubing 
just  downstream  of  the  venturi.  In  principle,  it  is  unchoked  under  any  other 
condition.  In  practice,  the  numerics  behave  badly  when  the  flow  depends  on 
both  upstream  and  downstream  pressure.  Smooth  calculations  must  either 
restrict  the  time  step  or  revamp  the  integration.  The  model  chooses  between 
two  venturi  flow  conditions,  fully  choked  or  no  flow.  The  effect  should  be 
small  in  the  upstream  chamber  which  is  the  target  of  this  analysis.  For 
designs  with  no  venturi  or  a  small  area  ratio  from  tube  to  venturi,  the  model 
is  inadequate. 


Burning  surface  comes  from  geometric  calculations  of  the  generator 
designers.  The  surface  is  essentially  constant  after  3  cm  burned  and  no 
longer  aominates  pressure  changes.  Expected  variations  in  area  after  the  3  cm 
should  affect  the  steady  pressure  level  only  by  a  factor  of  a  few  percent.  A 
coning  effect  may  increase  the  area  by  five  percent  which  translates  to  a  10? 
pressure  increase  IF  the  burning  rate  were  the  same  on  all  parts  of  the 
exposed  surface.  Slower  burning  at  the  edges,  due  to  heat  loss  through  the 
side  wall,  will  cause  the  coning.  The  net  result:  a  small  effect  on  the 
steady  pressure. 


Ignition  of  the  two  pellets  and  the  main  grain  need  not  be 
simultaneous.  A  reasonable  starting  condition  is  ignition  of  the  hot  pellet 
(AP  propellant)  and  both  chambers  filled  with  that  combustion  product  gas. 
Ignition  of  the  AN  propellant  (booster  pellet  and  main  grain)  can  then  be 
calculated  from  heat  transfer  to  the  surface  and  an  ignition  criterion.  Heat 
transfer  for  ignition  is  by  the  same  convection  and  unsteady  conduction.  The 
simplest  ignition  criterion  is  surface  temperature.  Since  there  seem  to  be  no 
ignition  data  for  the  propellant,  the  ignition  temperature  is  arbitrary. 
Temperatures  below  800K  assure  ignition  and  do  not  violate  widely-held 
theories.  At  higher  temperatures  ignition  depends  strongly  on  heat  transfer 
competition.  A  high  coefficient  helps  heat  the  propellant  faster  but  also 
cools  the  gas  faster.  The  race  frequently  goes  to  the  cooling. 


The  heat  transfer  coefficient  inside  the  combustion  chamber  is  itself 
uncertain.  Gas  velocity  is  low  near  the  grain  and  convective  heat  transfer  is 
inefficient.  But  the  grains  ignite  even  when  the  generator  fails  later.  The 
coefficient  must  then  be  high  enough  to  assure  ignition.  A  simpler  approach 
would  be  immediate  and  simultaneous  ignition  of  all  surfaces.  The  debate  is 
probably  academic  since  test  failures  seem  unconnected  to  ignition. 


The  uncertain  heat  transfer  affects  more  than  ignition.  Heat  transfer  to 
the  generator  head  and  filters  depends  intimately  on  the  coefficient.  The 
centrifugal  filter  will  have  a  higher  gas  velocity  and  thus  a  considerably 
higher  coefficient  than  the  gas  in  contact  with  the  propellant  surfaces.  The 
phenolic  filter  has  a  low  thermal  conductivity  and  therefore  heats  rapidly  at 
its  surface,  reducing  the  heat  transfer  quickly. 
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Data  obtained  from  various  sources  for  nominal  conditions  for  the 
generator  design  are  given  in  Table  1  (see  page  14).  Propellant  burning  rate 
measurements  showed  evidence  of  lot-to-lot  variation.  The  data  given  here  are 
for  a  ' new'  lot  . 


V.  ONE -DIMENSIONAL  FLOW  MODEL 

The  1  umped -parameter  model  above  is  based  on  assumptions  which  trade 
accuracy  for  simplicity  and  economy  .  However,  some  questions  about  system 
behavior  demand  a  more  detailed  analysis  .  Substantial  heat  loss  to  cold 
boundaries  can  lead  to  axial  property  variations  throughout  the  system  .  The 
behavior  of  both  the  nozzle  and  the  venturi  will  be  sensitive  to  the  local 
flow  properties  at  their  entrance  planes.  Since  the  system  failure  can 
apparently  be  reversed  by  addition  of  the  venturi,  it  is  important  to  explain 
how  this  device  alters  the  flow  field  . 

The  flow  field  is  assumed  one  -tiimens  ional  and  unsteady.  One -dimensional 
flow  is  a  reasonable  assumption  everywhere  except  in  the  combustion  chamber 
and  separator  chamber  (centrifugal  filter)  which  involve  low  speed  three- 
dimensional  flow.  To  avoid  this  complication,  the  flow  in  each  chamber  is 
node  led  as  a  reservoir  problem,  i.e.,  the  influx  is  assumed  to  stagnate  in  the 
chamber  and  the  outflux  is  accelerated  from  the  local  stagnation  condition. 
However,  a  further  complication  arises  when  the  usual  equations  of  motion  are 
applied  to  the  tubing,  venturi,  and  nozzle.  The  fact  that  both  the  venturi 
and  nozzle  have  large  values  of  entrance  area/throat  area  (33  for  venturi,  67 
for  the  nozzle)  requires  extremely  small  grid  spacing  within  these  devices  to 
resolve  the  flow  field,  particularly  if  the  operation  may  switch  between 
unchoked  and  choked.  Stability  then  forces  a  prohibitively  small  maximum  time 
step,  making  the  simulation  impractical.  A  reasonable  alternative  is 
construction  of  special  flow  field  elements  which  assume  inviscid,  quasi¬ 
steady  flow.  Note  here  that  quasi -steady  means  the  interior  flow  of  the 
element  responds  instantaneously  to  changes  at  its  boundaries;  it  does  not 
mean  steady-state  flow.  These  elements  are  explained  in  greater  detail 
below.  Finally,  combustion  products  from  both  kinds  of  solid  propellant  are 
treated  as  ideal  gases.  Mixture  properties  for  the  system  flow  field  are 
adjusted  (artificially)  in  proportion  to  the  mass  flow  rate  from  each 
propellant.  Ignoring  the  propagation  of  these  changes  along  streamlines. 

A.  Continuous  Flow  Field 


The  pipe  or  tubing  flow  field  is  described  by  the  solution  to  the  one¬ 
dimensional  unsteady  equations  of  motion  (accounting  for  wall  heat  loss)  given 
in  conservative  form  by 


t  +■  •  =  0 

t  z 


(15) 
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propellant.  If  there  is  an  age  effect,  the  margin  provided  by  the  venturi  may 
not  be  enough.  A  slightly  larger  venturi  (3%)  with  the  old  propellant 
failed.  The  safety  margin  cannot  be  estimated  from  the  data  and  analysis  to 
date.  On  the  plus  side,  the  propellant  has  long  been  used  in  other 
applications  with  no  strong  evidence  of  aging. 

4.  What  causes  the  pressure  difference  between  venturi  and  nozzle? 
Neither  model  calculates  it.  Nominal  model  input  calculates  unchoking  of  the 
venturi  just  after  peak  pressure  with  nearly  equal  pressures  thereafter;  but 
tests  typically  show  a  continuous  difference  of  about  0.7  MPa. 


VIII.  QUANDARY 

Is  the  probability  of  failure  of  the  present  design  high  enough  to 
justify  more  investment  in  understanding?  The  models  do  not  presently  qualify 
as  an  engineering  tool.  Continued  success  in  tests  will  make  the  problem  seem 
moot.  But  the  issue  of  age-related  change  will  not  be  answered  until  two 
years  pass  and  test  results  are  then  compared  to  early  production  results.  If 
and  when  the  first  test  failure  appears,  the  models  will  be  asked  to  steer  the 
engineering  of  the  repair.  But  only  if  the  manifest  inaccuracies  are  removed 
can  their  answers  be  trusted  to  be  any  better  than  the  crude  calculations  on 
which  the  design  already  rests. 


PRESSURE  (psi) 


TIME  (sec) 


Figure  11.  Pressure  Time-History  Predicted  by  One-Dimensional  Model. 

Incomplete  Combustion.  Solid  Line  =  Combustion  Chamber,  Dashed 
Line  =  Nozzle  Entrance  Plane. 

VII.  QUESTIONS  STILL  LURKING 

1.  How  did  heat  transfer  cause  the  failures?  Its  guilt  was  demonstrated 
in  tests  that  succeeded  when  the  valve  and  downstream  tubing  remained  at  the 
ambient  California  fall  temperatures  and  only  the  parts  upstream  of  the  valve 
were  conditioned  to  cold  temperature.  Simple  sensible  heat  loss  through  an 
ideal  equation  of  state  for  equilibrium  combustion  products  does  not  explain 
it.  Incomplete  combustion  may.  The  investigation  ended  without  a  useful 
post-mortem.  Also  unexplained  was  a  gradual  decline  in  pressure  for  about  25 
seconds  to  the  unacceptably  low  but  steady  pressure. 

2.  How  did  the  venturi  solve  the  problem?  The  basis  for  the  predictions 
of  the  venturi's  performance  is  attacked  by  the  finding  that  only  incomplete 
combustion  recreates  the  measured  pressure.  Venturi  design  calculations  were 
don-  with  equilibrium  composition  assumptions. 

3.  Will  the  fix  stay  fixed?  Propellant  aging  cannot  be  ruled  out  on  the 
evidence.  The  transition  from  success  to  failure  happened  with  two-year-old 
pro-ellant.  All  the  successful  tests  of  the  venturi  occurred  with  new 
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Figure  10.  Temperature  Time-History  Predicted  By  One-Dimensional  Model. 

Solid  Line  =  Exit  of  Dual  Chambers,  Dashed  Line  =  Nozzle 
Entrance  Plane. 


Because  heat  loss  seems  unable  to  explain  the  system  behavior,  attention 
focused  on  the  possible  lack  of  gas-phase  thermodynamic  equilibrium  as 
suggested  by  trials  with  the  lumped-parameter  model.  Setting  the  heat 
transfer  coefficients  at  their  nominal  values  and  then  imposing  a  fixed  state 
of  incomplete  combustion  (y-1  =  half  value, ~  three  times  nominal,  Tf  =  60/' 
nominal)  for  both  propellants  produces  the  results  shown  in  Figure  11  (also  on 
Figure  9  to  scale).  These  pressure  time-history  predictions  are  much  closer 
to  the  experimental  data,  although  the  simulation  again  insists  that  the 
venturi  will  unchoke  (at  t  -  0.45  sec)  after  maximum  chamber  pressure  is 
attained. 

The  dramatic  improvement  in  the  predictions  strongly  suggests  that  the 
gas-phase  combustion  products  may  be  undergoing  a  complex  "shifting" 
equilibrium,  possibly  with  condensation.  This  could  have  an  important 
influence  on  the  presence  of  a  shock  wave  in  the  venturi.  The  addition  of 
this  complicated  chemistry  to  the  flow  field  was  beyond  the  scope  of  the 
present  study. 
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GAS  TEMPERATURE  (Rankine) 


PRESSURE  (psi) 


system  behavior  is  the  same.  It  should  be  emphasized  that  doubling  the  heat 
loss  to  the  valve  and  tubing  which  separates  the  venturi  from  the  nozzle  is 
not  sufficient  to  keep  the  venturi  choked. 


Given  the  constraints  of  the  model,  the  simulations  demonstrate  that  only 
a  shock  wave  downstream  of  the  choked  venturi  is  capable  of  creating  the 
magnitude  of  total  pressure  loss  measured  in  the  actual  gas  generator. 
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Figure  9.  Pressure  Time-History  Predicted  by  One-Dimensional  Model.  Solid 
Line  =  Combustion  Chamber,  Dashed  Line  =  Nozzle  Entrance  Plane. 

However,  fluid  mechanics  coupled  with  the  assumption  of  gas-phase  equilibrium 
thermochemistry  predicts  the  venturi  will  not  remain  choked  after  maximum 
pressure  has  been  achieved. 

An  objective  fot  the  modeling  was  to  predict  the  influence  of  heat  loss 
on  the  flow  field  temperatures.  Figure  10  shows  a  comparison  of  temperature 
time-histories  for  the  two  cases  discussed  above.  The  solid  lines  represent 
gas  temperature  at  the  exit  of  the  dual  chambers,  and  the  dashed  lines  denote 
gas  temperature  at  the  nozzle  entrance  plane.  The  difference  between  these 
two  values  measures  the  effect  of  heat  loss.  Note  that  in  both  cases,  the 
temperatures  become  nearly  time  invariant  after  the  order  of  1  sec.  Thus, 
thermal  equilibrium  with  the  surroundings  is  established  much  sooner  than  the 
time  required  for  the  system  to  "fail." 
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PRESSURE  (MPa) 


Flow  equations  seem  an  unlikely  source  of  the  pressure  disagreement,  and 
burning  rate  was  well  characterized  by  independent  tests.  Heat  transfer 
offers  the  best  first  guess  for  the  difference  since  it  was  suspected  that 
heat  transfer  caused  the  problem  in  the  first  place.  But  reasonable 
variations  in  the  heat  transfer  descriptions  do  not  produce  a  credible  change 
in  the  peak  pressure  without  obliterating  some  other  aspect  of  the  problem. 
Enough  heat  transfer  to  reduce  the  pressure  peak  leads  to  gas  temperatures 
well  below  the  limits  of  equilibrium  thermodynamics. 

This  led  to  the  speculation  that  the  gas  thermodynamics  vary  from  the 
expected  equilibrium.  Some  incomplete  combustion  or  subsequent  condensation 
will  significantly  change  the  thermodynamic  properties  of  the  product  gas. 
Arbitrary  variations  in  thermodynamics  could  easily  be  made  to  test  hypotheses 
about  gas  chemistry  in  numerical  experiments.  Full  and  defensible  chemical 
effects  are  left  for  another  day.  Physical  evidence  for  a  non-equilibrium 
condition  comes  from  window-bomb  tests  by  the  Jet  Propulsion  Laboratory  (Leon 
Strand).  Films  of  the  combustion  show  hot  ash,  a  less  visible  flame,  and  i 
larger  ash  residue  at  lower  pressures. 

For  thermodynamic  consistency,  a  mixture  of  gases  was  assumed  to  consist 
of  fully  reacted  equilibrium  products  and  an  arbitrary  intermediate  product  of 
incomplete  combustion.  The  progress  variable  that  controlled  the  conversion 
of  intermediate  to  final  products  could  also  be  arbitrary.  Only  simple  linear 
time  dependence  was  tried.  Trial  and  error  variations  in  y,  n  ,  Tf  were  made 
until  calculated  pressures  matched  test  results.  The  values  to  produce  this 
agreement  were:  (y-1.)  half  nominal  value, ^  three  times  nominal,  and  T^  60? 
of  nominal.  Time  dependence  of  the  progress  variables  ended  arbitrarily  at  2 
sec . 


Ignition  of  only  part  of  the  propellant  surface  is  another  candidate  for 
error  source.  If  the  burn  rate  exponent  (Eq.  (10))  is  0.45,  a  pressure  drop 
by  a  factor  of  three  would  need  surface  area  of  roughly  half  the  geometric 
area.  Given  the  highly  gaseous  igniter  pellet  products,  a  half  ignited 
surface  seems  unlikely. 

B.  One-Dimensional  Flow  Model 


Figure  9  shows  a  comparison  of  pressure  time-histories  predicted  by  the 
one-dimensional  model  for  three  different  cases.  The  solid  lines  represent 
combustion  chamber  pressure,  while  the  dashed  lines  denote  static  pressure  it 
the  entrance  plane  to  the  nozzle.  Two  important  results  are  demonstrated  bv 
the  predictions  for  the  nominal  case.  First,  maximum  chamber  pressure  is 
greater  than  7  Kpsi,  in  the  range  predicted  by  the  lumped-parameter  model. 
Second,  the  venturi  unchokes  soon  after  the  combustion  chamber  achieves 
maximum  pressure,  e.g.,  approximately  0.2  sec  in  this  case.  (When  tie  dn-.be  : 
and  solid  lines  become  coincident  on  the  scale  of  the  plot,  the  venturi  L 
operating  subsonic.)  Roth  features  are  in  direct  conflict  with  the 
experimental  data.  Because  of  the  uncertainty  associated  with  some  .isn-u  *s 
the  heat  transfer,  the  case  was  rerun  after  setting  all  heat  transfer 
coefficients  to  twice  their  nominal  values.  These  results  (labeled  ,i. i ,, 
Heat  Loss)  are  shown  in  Figure  9  and  confirm  the  conclusion  from  the  lumped 
parameter  model  that  system  behavior  is  relatively  insensitive  to  t  ha  ores  ■ 
heat  transfer.  The  venturi  remains  choked  for  a  longer  time,  'mi*  the  ,>en,  r.n  ' 


V  ■ '.7*  ““ “ 


in  the  model  results.  Choking  in  the  venturi  is  one  possible  explanation  as 
is  a  pressure  flow  loss  in  the  tubing.  Unfortunately,  neither  explanation  is 
supported  by  model  findings  or  by  fluid  mechanics  intuition. 


Figure  7.  Typical  Cold-Temperature  Test  Pressures 

At  long  times  the  predicted  and  measured  pressures  agree.  Figure  8  shows 
calculated  pressures  out  to  65  seconds  for  the  upstream  chamber. 


2 


[»T»Tt, 


Q-  4000 


3  3000 

iu  // 

o-  2000  A  < 


* 


1000 


0  0.2 


0.4  0.6  0.8 

TIME  (s) 


Figure  5.  Predicted  Pressures  at  Early  Time 


2000i 


h+-T 


1000 


VENTURI  ENTRANCE 


NOZZLE 


S  500tt 


0  0.5  1.0 

TIME  (s) 

Figure  6.  Test  Pressures  at  Early  Time 


Figure  7  shows  a  typical  cold-temperature  test  history.  The  results  also 
miss  the  measured  difference  between  venturi  entrance  and  nozzle  pressures.  A 
consistent  difference  gap  of  at  least  0.7  MPa  appears  in  the  test  data  but  not 
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The  jump  equations  for  flow  through  the  area  change,  3+4,  are  similar  to  the 
above  set.  The  system  is  completed  by  the  addition  of  Eq.  (19)  written  along 
the  left-running  characteristic  line  which  intersects  point  4;  this  provides 
the  communication  link  to  the  flow  in  the  downstream  tube  and  nozzle.  The 
total  system  of  equations  was  hand  reduced  to  four  equations  in  four  unknowns 
and  solved  with  the  stiff-equation  root  finder. 

The  equations  governing  the  unchoked  (subsonic)  operation  of  this  element 
are  similar  to  the  above,  without  the  complication  of  the  normal  shock  wave. 

The  crucial  decision  whether  the  venturi  element  is  choked  or  unchoked  is 
based  on  the  solution  to  he  complete  equation  system.  If,  during  the  choked 
solution,  the  upstream  shock  Mach  number  is  driven  below  unity,  the  solution 
procedure  switches  to  the  unchoked  equations  and  retries  the  solution.  If, 
during  the  unchoked  solution,  the  effective  choked  area  [implicit  in  Eq .  (31)] 
falls  below  the  geometric  throat  area,  a  switch  is  made  to  the  choked 
equat  ions . 

Convective  heat  transfer  to  the  steel  valve  and  tubing  assumes  fully- 
developed  turbulent  pipe  flow;  the  heat  transfer  coefficient  [for  Eq.  (5)]  is 
given  hv 


h  =  0.025 
c 


Pr°*4  (ReD)0,8 


The  massive  steel  valve  is  assumed  to  be  an  infinite  sink,  but  the  outer 
surface  of  the  tubing  is  allowed  to  radiate  as  a  black  body  to  the  cold 
ambient  temperature.  A  compromise  solution  in  the  combustion  chamber  assumes 
heat  loss  to  infinitely  thick  steel  walls  at  a  rate  ten  times  that  for 
stagnant  flow,  i.e.,  Nusselt  number  =  20. 


VI.  RESULTS 

A.  Lumped- Parameter  Model 

For  the  nominal  design  and  conditions  (see  Table  i)  the  predicted  earlv 
pressure  is  shown  in  Figure  5,  Peak  pressure  far  exceeds  the  measured  peak  as 
seen  in  Figure  b,  although  the  time  of  the  peak  nearly  coincides  with  test 
lata.  Time  to  peak  depends  weakly  on  pressure  because  it  depends  on  burning 
rate  of  the  pellet  which  itself  has  a  low  pressure  dependence. 


Figure  4.  Schematic  and  Notation  for  Venturi-Valve  Element 

For  the  case  when  the  venturi  is  choked,  the  solution  for  the  entrance 
flow  (1)  is  uncoupled  from  the  flow  downstream  of  the  throat  and  can  be  found 
in  the  manner  described  in  paragraph  C  above  for  the  choked  nozzle.  In  the 
region  downstream  of  the  throat,  Eq.  (28)  becomes  a  relationship  between  the 
area  location  of  the  normal  shock  wave,  A„w,  and  the  upstream  Mach  number, 
Msl,  ie., 

Y  +  1 

A  =  {A  /C  M  }  [1  +  X-=-i  M  2]2^Y  "  l)  .  (30) 

sw  tv  f  SI  2  SI 


Further  algebraic  manipulation  provides  a  transcendental  expression  between 
the  downstream  shock  Mach  number,  M^,  and  the  exit  plane  Mach  number,  M2 
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(31) 


When  the  Mach  numbers  are  determined,  total  pressure  loss  across  the  shock 
wave  and  the  new  stagnation  conditions  are  simple  to  compute.  The  equations 
describing  flow  in  the  valve,  accounting  for  heat  loss  to  the  boundary,  are 
given  by 
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Figure  3.  Schematic  and  Notation  for  Choked  Nozzle  Element 


Tsentropic  flow  provides  a  unique  relationship  between  area  ratio  and  Mach 
number,  viz. 

Y  +  1 

O  ( y  »  1  } 

choked  area  _  A_  _  ^  M  jj  +  Y  -  1  m2j  (28) 

local  area  A  f  2 


where  -  ( 
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Y  +  1 
2 ( Y  -  1) 


Given  the  .area  ratio,  Eq.  (28)  is  a  transcendental  expression  for  M  which 
chances  only  when  gas  composition,  Y,  changes.  At  the  entrance  plane,  Eq. 
(28)  determines  Mj  which  also  must  satisfy 


o  1  u i 2  -  YP^2  =  0  .  (29) 

This  relationship,  along  with  Eqs.  (18)  and  (20)  evaluated  at  the  entrance 
plane,  uniquely  determine  the  flow  properties  at  1.  It  is  important  to  note 
that  this  choked  nozzle  solution  responds  (instantaneously)  to  any  changes  in 
entrance-plane  properties;  it  does  not  enforce  a  time-independent  value  of 
choked  mass  flow  rate. 

D.  Venturi -Valve  Element 


Construction  of  this  element  follows  that  of  the  nozzle  but  is  more 
•or.pl  ex.  It  consists  if  a  converging  and  diverging  section  which  is  attached 
t  i  a  constant-area  valve  terminated  by  an  area  change  to  match  the  diameter  o 
the  'ubing  section  leading  to  the  nozzle  (see  Figure  4).  Depending  upon  the 
over  ill  static  pressure  difference,  this  element  can  operate  unchoked  with 
subsonic  flow  throughout,  or  choked  with  a  normal  shock  wave  standing 
s  'nowhere  in  the  divergent  section.  In  both  cases,  the  flow  field  is  assumed 
quas i -s t eadv  and  isentropic,  with  provisions  for  the  shock  wave  (if  present) 
ami  wall  heat  loss  in  the  valve. 
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qwc  =  heat  loss  through  wall  area  Awc  of  combustion  chamber 

Obviously,  chamber  pressure,  Pc,  follows  from  the  equation  of  state,  P=P(p,e). 

The  separator  chamber  contains  no  combustible  material  and  hence  has 
constant  volume,  Vgc.  The  remaining  equations  are: 

mass 
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For  low  speed  flow,  the  mass  flux  between  chambers  can  be  written  as 

[e. g. ,  see  Ref.  2,  p.  95  ] 
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RT 
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where  6  is  the  solution  to 


3  (1  ♦  3/4) 


(27) 


with  a  similar  construction  for  Wjj-  The  influence  of  the  pipe/tubing  flow  is 
communicated  by  the  compatibility  condition  [Eq.  (19)]  along  the  left-running 
line  reaching  the  plane  T2.  With  this  addition,  the  dual-chamber  flow  is 
uniquely  determined.  "Tie  resulting  equation  system,  however,  is  quite  stiff; 
a  special  stiff-equation  root  finder  is  required  for  solution. 

C.  Choked  Nozzle  Element 


This  special  solution  element  (see  Figure  3)  assumes  that  the  nozzle 
remains  choked,  the  flow  field  is  quasi-steady,  and  isentropic.  The 
implication  of  quasi-steady  is  that  the  ratio  of  nozzle  length  to  local  sound 
speed  is  much  less  than  the  time  required  for  a  change  in  the  system  flow 
field.  This  should  be  a  good  assumption  here. 


"A.H.  Shapiro,  Compress ible  Fluid  Flow,  Vol.  I,  Ronald  Press,  1953. 
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Figure  2.  Schematic  and  Notation  for  Combustion  Chamber  Model 


Low  speed  flow  (M^<<1)  is  a  reasonable  assumption,  as  are  spatially 
uniform  pressure  and  zero  mass-average  velocity  in  each  chamber.  Real-world 
flow  losses  in  transfering  mass  between  chambers  and  into  the  tubing  entrance 
are  modeled  with  "orifice-loss"  coefficients. 


For  the  combustion  chamber: 
volume  change 
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where  'Hq  j  =  mass  flux  through  orifice  area 

hT',  hTT°  =  flame  enthalpies  of  the  two  propellants 
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momentum 
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where  w  =  pu, 


£=  P  (e  +  u2/2) 


(17) 


and  P=P(p,e)  is  prescribed  by  the  equation  of  state.  The  numerical  solution 
is  predicted  with  MacCormack's  explicit  scheme.  Coupling  between  this 
unsteady  solution  and  the  special  flow  field  elements  discussed  below  is 
accomplished  with  method-of-character ist  ics  compatibility  conditions  along 
characteristic  directions,  all  of  which  follow  from  the  above  equation  system 
written  in  characteristic  form: 

dP  +  pa  du  =  -Qdt  (18) 
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B .  Combustion  Chamber  and  Separator  Chamber 

The  dual  chamber  problem  is  modeled  as  two  reservoirs  connected  bv  an 
orifice,  Cl.  The  tubing  inlet  is  an  outflux  boundarv,  T2,  to  the  separator 
chamber  (see  Figure  2). 
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local  sound  velocity  [see  Eq.  20] 
cross  sectional  area 

surface  area  of  main  propellant  grain 
nozzle  throat  area 
venturi  throat  area 
propellant  surface  area 

specific  heat  at  constant  pressure,  volume 

orifice  loss  coefficient  [see  Eq.  26] 

diameter  of  circular  cross  section  (2r„) 

w 

specific  internal  energy 

total  energy  per  unit  volume 

defined  in  Eq.  28 

specific  enthalpy  (e  +  P/p) 

reference  enthalpy  (heat  of  formation) 

convective  heat  transfer  coefficient 

thermal  conductivity 

axial  length  of  valve 

Mach  number 

molecular  weight 

static  pressure 

Prandtl  number 

boundary  heat  flux 

solid  propellant  regression  rate 

wall  radius  of  circular  cross  section 

Reynolds  number  based  on  diameter 
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t 


t  ime 


T  temperature 

u  velocity 

v  volume 

z  axial  distance 

Zv  volume  burned,  defined  in  Eqs.  21a  &  21b 

y  isentropic  index 

P  density 

oj  mass  flux  (pu) 

(  )t  partial  derivative  wrt  time 

(  )  partial  derivative  wrt  distance 

(  )  pertaining  to  combustion  chamber 

(  )RC  pertaining  to  separator  chamber 

(  )  pertaining  to  valve 

(  )j  pertaining  to  propellant  1  (Ignition  Pellet) 

II  (Main  Grain) 
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